Control of Josephson current by Aharonov-Casher Phase in a Rashba Ring 
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We study the interference effect induced by the Aharonov-Casher phase on the Josephson current 
through a semiconducting ring attached to superconducting leads. Using a ID model that incor- 
porates spin-orbit coupling in the semiconducting ring, we calculate the Andreev levels analytically 
and numerically, and predict oscillations of the Josephson current due to the AC phase. This result 
is valid from the point contact limit to the long channel length limit, as defined by the ratio of the 
junction length and the BCS healing length. We show in the long channel length limit that the 
impurity scattering has no effect on the oscillation of the Josephson current, in contrast to the case 
of conductivity oscillations in a spin-orbit coupled ring system attached to normal leads where im- 
purity scattering reduces the amplitude of oscillations. Our results suggest a new scheme to measure 
the AC phase with, in principle, higher sensitivity. In addition, this effect allows for control of the 
Josephson current through the gate voltage tuned AC phase. 

PACS numbers: 73.23.-b, 72.10.-d, 74.45. +c, 74.50. +r, 03.65.Vf 



I. INTRODUCTION 

When a quantum system undergoes cyclic evolution 
motion in parameter space, its wave function acquires 
a geometric phase which strongly influences the trans- 
port property of the system. The best known example 
of such phase is the Aharonov-Bohm (AB) phase which 
is the relative phase shift between two charged particle 
paths enclosing a magnetic flux. In the last few years, the 
Aharonov-Casher (AC) phase, 1 another example of a ge- 
ometric phase, has been studied both experimentally 2,3 
and theoretically^^ in spin-orbit (SO) coupled meso- 
scopic semiconducting rings. As the charge-spin dual 
of the AB effect, the AC phase in semiconductors has 
attracted great attention because it can be easily con- 
trolled by a gate voltage. So far, most of the litera- 
ture has focused on Rashba SO coupled rings attached to 
normal metal leads. The AC phase has been confirmed 
in these type of systems by observing the conductance 
oscillation as a function of the gate voltage that is re- 
lated to the SO interaction (SOI) strength. However, 
due to the large background current, the typical ampli- 
tude of the oscillation is no more than 10% of the total 
observed conductance^ The AC phase has also been dis- 
cussed for magnetic vortices in type-II superconductors 
and Josephson junction array i 7 ^^^ Recent studies in the 
effect of SOI on Josephson curren t 11 i 12 ' 13 i 14 ' 15 has focused 
on heterostructures comprised of a two-dimensional elec- 
tron gas (2DEG) or a quantum dot placed between 
two superconducting leads. A natural and interesting 
step is to combine the AC and Josephson effects to 
explore new regimes in semiconductor-superconductor 
structures, which has not been explored to our knowl- 
edge. 

In this article, we study the AC effect in a Joseph- 
son current traversing a ring with Rashba type SO in- 
teraction. The Josephson current in a superconducting- 
normal-superconducting (S/N/S) junction mainly origi- 



nates from the Andreev reflect ion t 16 i 17 (electrons being 
retroreflected to holes) in two Superconducing/Normal 
(S/N) interfaces. The excitation spectrum below the 
superconducting gap consists of discrete energy levels, 
called Andreev levels J£ Andreev levels are affected by 
both the phase and the transmission coefficient of elec- 
trons and holes in the normal regionj 19 ' 20 i 21 We study 
here the effects of SOI on the Andreev levels and the 
corresponding Josephson current through the ring, illus- 
trated in Fig. [TJ both analytically and numerically. Our 
calculations predict oscillations of the Josephson current 
controlled by the AC phase. This effect can be used as 
an alternative sensitive way to detect the AC phase. Be- 
cause the dephasing electrons will not contribute to the 
current in a Josephson junction but will contribute to 
the current in a normal junction attached to two non- 
supercoducting leads, the observed amplitude of the cur- 
rent oscillations due to AC effect in a ring Josephson 
junction should be much larger than the conductance os- 
cillations in the normal ring junction. 




FIG. 1: Two superconductor leads are attached at the points 
ij) = and ip = 7T separately on the ring, Josephson current 
can be obtained in this junction. 

The outline of this paper is as follows. In Sec. II, we 
introduce our theoretical model of the Josephson current 
in a 1-D mesoscopic SO coupled ring-shaped Josephson 
junction. In Sec. Ill, we use the boundary conditions 
of a multiple terminal junction and present the numeri- 
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cal calculations of the Andreev levels and the Josephson 
currents in limits not explored by the analytical solu- 
tion. In Sec. IV we present an experimental set-up for 
the observation of the effect and in Sec. V summarize 
our conclusions. In Appendix A, we present the details 
of the sub-gap states of quasiparticles in superconduc- 
tor. In appendix B we show the details of calculating the 
eigenstates and the eigenenergies of the ring in the tight 
binding model. 



II. ANALYTICAL DISCUSSION OF ANDREEV 
SPECTRUM DUE TO SPIN-ORBIT COUPLING 

The Andreev bound state is a kind of electron-hole 
bound state in the normal region of the Josephson junc- 
tion depicted by Fig. [2] In the superconductor leads, the 
excitation spectrum consists of the positive eigenvalues 
of the Bogoliubov equation 
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where H(r) is the one-electron Hamiltonian, /i is the elec- 
trochemical potential, U(r) is the scalar potential, A is 
the superconducting pair potential and m* is the effective 
mass. When the excitation energy is less than the gap 
|A|, the quasiparticles of Eq. (TT|) in the normal region 
will be reflected by the pairing potential and form the 
Andreev bound state. Unlike the bound states in a usual 
square well, the Andreev levels carry electric current 
which contributes to most of the Josephson current.— 
Generally speaking, the bound states can be described in 
such a way that if we have a right moving subgap parti- 
cle on the left S/N interface Fig. [21 after gaining a phase 
described by the matrix P ir due to the particle propa- 
gation from left to right, reflection on the right interface 
described by the scattering matrix S r , a further phase de- 
scribed by the matrix P r \ from the motion from right to 
left, and a reflection on the left interface with the scatter- 
ing matrix Si, the state should come back to its original 
state. This is true only when Det(l — P r iS r Pi r Si) = 0, 
which determines the quantized Andreev levels. 

In zero magnetic field and zero SOI, neglecting Fermi 
wavelength mismatch and the barrier on the S/N inter- 
face and assuming that the transmission function of elec- 
trons in normal region is energy independent, these dis- 
crete energies are determined by^ 
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where £ = Hvp/2A is the BCS healing length, <fr 
4>2 — 4>\ is the superconducting phase differences, n 
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FIG. 2: A step- like superconducting pair potential is assumed 
in the Josephson junction. The right moving particle in the 
normal region will come back to its initial state after two 
Andreev reflections. The change of the color of the arrow 
after each reflection means an electron is reflected to a hole 
and a hole is reflected to an electron. 



0, ±1,±2---, 9 is the additional phase due to a point 
impurity potential VS(x) in the normal region, a is the 
distance between the point impurity and left interface, 
and T and R are corresponding to the transmission prob- 
ability and reflection probability due to a point impurity 
when energies are close to the fermi energy. The term 
(EL/ A£o) is approximately the phase shift acquired from 
free electrons and holes propagation in the normal re- 
gion. Eq. ([3]) demonstrates that the Andreev level will 
be determined by both the phase shift and the transmis- 
sion function of electrons and holes in the normal region. 
An important limit is when the length of the junction 
is much shorter than the BCS healing length (L <C £0) 
so that 7 ~ 0. In this case, the Andreev level takes the 
formal 



£ = Acos(-)=Ayi-Tsin 2 (|). 



(4) 



Next we focus on the Josephson current in a one- 
dimentional (1-D) Rashba SOI ring; this current can be 
computed from the superconducting phase dependence 
of the Andreev bound states. For our model calculations 
we assume a step-like superconducting pair potential A, 
which is zero in the Rashba ring: 



A e-^/ 2 , 
{0, 
A e^ 2 , 



left of the ring 
in the ring 
right of the ring 



(5) 



where +(— )<^/2 are the superconducting phase in the 
right(left) leads. 

Because the SOI keeps time-reversal symmetry, the co- 
efficient of the Andreev reflection is almost the same as 
without SOI. As a result, Eq. ([3]) is still valid qualita- 
tively for a ring-like junction with SOI. The ID Rashba 
Hamiltonian for a 1-channel ring of radius ro without 
magnetic field reads ^ 
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where r$ is the radius of the ring, ljq = h/{m*r^) with m* 
being the effective mass of the electron in the ring, and 
uor = 2a/ (hr ) with a being the the strength of the spin- 
orbit coupling. The cigcnfunctions of this Hamiltonian 
take the form: 



= exp(m^) ( _ s "° n S ( ^| e , 
The associated eigenenergies read 
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where n e is the quantized angular momentum number of 
the electron with spin along —z direction, m e is the quan- 
tized angular momentum number of the electron with 
spin along z direction, /i is the chemical potential of the 
system and C e — 2 (^ + ^) is a parameter independent of 
SOI. 

Since the hole is the time-reversal state of the electron 
and SOI keeps time-reversal symmetry, we have = 
—H e . Therefore, the hole states are the same as the 
electron states and the associated eigenenergies take the 
similar form: 
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We next consider the phase difference between to the 
electron and hole propagation in the normal region. The 
phase difference in the upper part of the ring for the 
spin-up and spin-down quasiparticles after being Andreev 
reflected two times take the form 



In = \n e - n h )TT = 2 — ^ 7T, 



7m = (m e - mh)7r 



and are therefore identical. Here 7„( m ) is the phase shift 
of the quasiparticle with spin along — z(+z) direction. 
This identity reflects the fact that there is no spin split- 
ting of the phase shift in this case and is consistent with 
the time-reversal symmetry of SOI. The phase acquired 
in the Andreev reflection from the interfaces lead to a 
different phase shift since the spin up electrons and spin 
down holes have different momentaJ^ However this phase 
shift splitting is very small and will not affect the zeros 
of the Josephson current. Therefore we will not consider 
this splitting in our discussion. 

The effect of the Rashba interaction on the trans- 
mission function has been studied previously with a 
conductance oscillation due to the AC phase confirmed 
experimentally 2,3 and theoretically^^ Under the as- 
sumption of a perfect coupling between the leads and 
the ring and neglecting backscattering, the transmission 
function takes the formi 
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When considering the short junction limit, L <C £c>! 
the energy E of the Andreev bound state is affected by 
the normal-state transmission T through Eq. J4{. The 
Josephson current in the low temperature limit takes the 
form 21 
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Because the transmission probability is affected by the 
AC phase, an oscillation of the Josephson current due to 
SOI should also be observed. This oscillation will be dif- 
ferent from the conductance oscillation because even in 
the short junction limit, the current is nonlinearly depen- 
dent on the transmission function. Although this conclu- 
sion is obtained in the short junction limit [L <gC £ )j since 
the zeros of the transmission probability are only depen- 
dent on the SOI, a similar oscillatin can be expected at 
any value of £/£o- We will show this to be the case in 
our numerical calculations. 



III. NUMERICAL CALCULATION OF 
ANDREEV LEVEL AND JOSEPHSON CURRENT 



In this section, we present our numerical investiga- 
tion on the relation of the Josephson current to the SOI. 
The effetive tight-binding model Hamiltonian of the ID 
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Rashba ring is given as 5 *: 



A' 
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where ipj = 2n(j - l)/N, <p jd+1 = (<pj + <fj+i)/2, 
Sip = 2ir/N, t so — a/2a, r is the radius of the ring and 
to = ft 2 /2ma 2 with a being the lattice constant of the 
tight-binding model. The length scale a is not related 
to the lattice constant of the material but is an artifical 
length scale used in modeling the continum Hamiltonian, 
Eq. ([6]) , and all physical results obtained should be inde- 
pendent of this length scale; e.g. the Fermi energy scales 
considered should be near the bottom of the tight-binding 
Hamiltonian band, Ep <C 4to- The eigenf unctions of this 
Hamiltonian are the same as in Eq. 10 . The eigenenergy 
is obtained as 



E 

Eem 
Ehn 

where 
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n e (h) = (±K(h) ~ 0)/5<P + I/ 2 , 
m e(h) = (±A e (/j) + (3)/5ip- 1/2, 
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cos{8ip /2) 



=). The detailed derivation 



Here 3 = arccos(- , 

of the result shown above is presented in Appendix B. 
Defining the AC phase as (f>Ac = w so tt = (2/3/Sip — l) and 
comparing n and m in the tight binding model with those 
in the continuous case, we can find that the term w so tt = 
(2(3/5<p — 1) is the counterpart of the term (a/1 + Q 2 R — 
1)tt in the continuous case and they can be shown to be 
equal in the limit a — > 0. In the following discussion, we 
will often use the dimcnsionless parameter w so instead of 
tso/t in our numerical calculations and plots as a measure 
of the SOI strength since the AC phase is a monotomic 
increasing function of it. 

If we ignore the disoreder in the ring, the phase shift 
matrices Pi r and P r i are diagonal matrices and each diag- 
onal element takes the form exp(i(n e (/j) — n^e))^)- The S 
matrices 5; and S r are calculated by solving the bound- 
ary condition at each joint of the superconductor leads 
and the ring. In a ring system, the joint of the lead 
and the ring can be considered as a 3-way junction. The 
boundary condition of this 3-way junction in the tight- 
binding model is obtained from the principle that the 



wave function at the joints is also the eigenfunction of the 
system. 25 This boundary condition in the tight-binding 
model allows us to consider more realistic cases such as 
having different effective mass of electrons in the super- 
conducting leads and normal region. Once S matrices 
Si and S r are obtained, and substituting them into the 
bound state condition— 



Det(l - PrlSrPwSl) = 0, 



(18) 



we can find the relation among the Andreev level E, su- 
perconductor phase difference </> and the SOI strength. 
Eq. (fTS)) contains all of the Feyman paths and allows us 
to consider the effect of multiple backscattering within 
the ring. 

In the previous section, to get the analytical result, 
we made the assumption that the coupling between the 
ring and the leads are perfect and backscattering was 
therefore neglected. However, this idealization is not 
a good approximation in general since there is usu- 
ally a very large Schottky barrier on the superconduc- 
tor/semiconductor interface and the backscattering can 
also strongly affects the electron transmission function 
through a ring. All of these factors are considered in our 
numerical calculation. 

First, let us focus on the effects of the backscatter- 
ing. We consider the limit L <C Co- I n the tight bind- 
ing model, Co = hvf /2A — ato sin(fc/a)/A. We assume 
A = O.OOlio, Ef — to (i.e. one quarter of the band- 
width), Co = 500\/3a and the perimeter of the ring is 
100a which is smaller than the healing coherence length 
Co- We choose Ef = to because the transmission proba- 
bility of the electron through the ring oscillates slower in 
this energy range compared to its zero minimal value and 
when the energy is close to the bottom of the band. As 
a result, the electron's behavior in this energy range sat- 
isfies the approximation we made in the previous section 
and more directly confirms our analytical discussion. The 
parameters corresponding to the experimental system 2 
will be discussed at the end of this section. 

The Andreev levels are presented in Fig. [3] as a func- 
tion of the superconducting phase difference <fi f° r dif- 
ferent SOI strength in the dimensionless variables E/A, 
w so , and 4>/ir. When the AC phase is equal to (2n+ 
the energy of the Andreev states is A and independent 
of (f>. Since the Josephson current is proportional to the 
first derivative of the energy E of the Andreev bound 
state with respect to the phase difference cf>, this implies 
that the zeros of the critical current occurs at the values 
(2n + l)ir of the AC phase. This conclusion is consis- 
tent with our theoretical result in Eq. ((3]). In Fig. 21 we 
show the Andreev levels and the normalized Josephson 
currents, (l/Ad)E/d<fi — —I / (2eA/H) 1 associated with 
the different AC phases vs. superconducting phase dif- 
ference. The results show that the Josephson currents 
amplitudes are tuned by the AC phase. Under the short 
junction limit, the energy of the Andreev levels is always 
equal to the superconducting gap A when the phase dif- 
ference = and minimum when (f> — n for any SOI 
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FIG. 3: Top: Andreev levels vs. superconducting phase 
difference 4> for different SOI strength measured by AC 
phase(ui so 7r). Bottom: Andreev levels (color plot) vs. su- 
perconducting phase difference, <f>, and SOI measured by w so - 
Here N = 100, E s = to, A = 0.001t , and £ = 500 V3a. 



strength. Therefore the value of E(0) — E(ir) is a di- 
rect indicator of the amplitude of the critical current. 
E(ir) — E(0) versus AC phase, 4>ac — w so ir, is shown in 
Fig. O The zeros correspond to the value (2n + 1)tt of 
the AC phase and the peaks are not always correspond- 
ing to 2n-K but show a doubled peaked structure centered 
around 2mr. This is because of weak localization effect^ 
similar with the AC phase effect on conductance oscilla- 
tion in normal junctions.—^ 

We next consider the long channel length limit, L > £o, 
where not only the transmission function but also the 
phase shift due to the AC effect will affect the Andreev 
levels. As an example, Fig. [6] plots the Andreev levels 
as a function of superconducting phase difference <f> with 
different AC phases in a long-junction. Here the number 
of sites in the ring is changed from 100 to 2200, i.e. the 
length of the ring changes by a factor of 11, while the 
other parameters are unchanged. The Andreev levels are 
shifted due to SOI and the Josephson current vanishes 
when w = (2n + l)w as in the case of the short-junction 
and confirms that the backscattering will not affect the 
zeros of the critical current in the Josephson junction in 
both the short and long channel length limits. 

When the Schottky barrier and the Fermi wavelength 
mismatch are also considered, the current-SOI relation 
will be more complicated because the electon transmis- 
sion function will be strongly affected by these factors. 
However our numerical resutls show that the zeros of 
the critical current are the same to the analytical result 
Eq. (fT"4|) . We explore the Josephson current oscillations 
due to the AC phase at the lower, more realistic, Fermi 
energy E f = 0.1t , and take A = 0.0002i , Co = 500^0. 
The density of the electron in the experiment using HgTe 
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FIG. 4: The Andreev levels (left), E/A, and normalized 
Josephson currents (right), (l/Ad)E/d<f>, vs. supercon- 
ducting phase for different SOI strength measured by AC 
phases(w so 7r). Here N = 100, E f = to, A = 0.001t , and 
Co = 500\/3a. 
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FIG. 5: Difference of highest and lowest Andreev levels vs. 
SOI strength measured by w so . The difference is proportional 
to the critical current. Here N = 100, Ef = to, A = O.OOlto, 
and Co = 500\/3a. 



hetero-structures varied from — 1-85 x 10 cm to 
ri2D = 2.21 x 10 12 cm~ 2 and the effective mass is m* — 
Q.Oilmo^ To make our parameters match the experimen- 
tal data, by using the relation Et = O.Ito = h 2 /2m*a 2 = 
h 2 k 2 /2m* , Trri2D — k 2 , the parameter a is required to be 
around 1.5 nm. When we choose N = 2200, the perime- 
ter of the ring in our calculation is 3.3/im which is the 
same order to the ro = l//m in the experiment.— The su- 
perconducting gap A = 0.0002i = 0.002S/ » 1.72 K 
which is the same order of the gap of the conventional 
superconductor such as Al where T c = 1.18 K.— To cap- 
ture the effects of the Schottky barrier and the Fermi 
wavelength mismatch on the short and long Josephson 
junction more clear, we explore three ring length scales 
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FIG. 6: Andreev levels as a function of superconducting 
phase difference (f> for different SOI strength measured by AC 
phase(«) so 7r). Here N = 2200, Ef = to, A = 0.001t , and 
Co = 500\/3a. 
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L = 200a, L = 2200a, and L = 10000a, shown in 
Fig. [7J The Schottky barrier is chosen to be 0.8to and 
the Fermi wave vector mismatch is kt/k s f — 0.2 where 
kf and k s f are corresponding to the Fermi wavevector 
in the normal region and superconducting leads. The 
above three oscillations have a common characteristic, 
when the AC phase is equal to (2n + l)7r, the critical 
current is zero. This means the AC effect is present even 
at a large Shottky barrier and Fermi wavelength mis- 
match. In the bottom panel of Fig. [71 the critical current 
has negative values, which would indicate that in this 
case there could exists a ir junction^ This phenomena 
is due to the large difference of the phase shifts acquired 
by the electron and hole in the normal region which 
is related to the nonzero momentum Q of the cooper 
pair o 36 ' 38 and has been observed in the superconduct- 
ing/feromagnetic/supconducting junctions i 39 ' 40 ^ 41 How- 
ever the 7T junction can not be realized in our model be- 
cause it requires the condition L»^o and the electron 
will be totally decohered and the Josephson current will 
not survive. 



FIG. 7: Normalized criticla current vs. the SOI strength 
measured by w ao in different lengtgh scale of the SOI ring. 
Here the Fermi energy Ef = O.lto, the Fermi wavelength mis- 
match between the SOI ring and the superconducting leads is 



P- = 0.2 and N = 200. 



given byj^ 

H = H r + H l + Ht 

H T = ^(TkqlC^Ckl +T kql C^C k l + T kq1 C k1 C Ql 
kg 

+n q iC^c qL ) (19) 

where Ck(C^) is the annihilation(creation) operator in 
left side, C g (C+) is the one in right side, and T^q] is 
the tunneling amplitude of the spin up electron from the 
left lead with the momentum k to the right lead with 
momentum q. In the short junction limit L <C £, we 
can assume T^ q to be a constant T, and the Josephson 
current is calculated frorn^ 3 . 



The advantage of observing the AC effect in the ring 
Josephson junction lays on the fact that the dephasing 
electrons do not contribute to the Josephson current. 
The dephasing effect due to the impurities can be consid- 
ered through the Green's function method. The tunnel- 
ing Hamiltonian used in calculating Josephosn current is 



1 = — / diolm 
h 



(T£( W )T*&H) U } 



(20) 



where D r {u) = Det[l- \T\ 2 T 3 g r L (u)T 3 g r R (uj)] and r 3 is the 
usual Pauli matrix. This method gives the same analyti- 
cal form of the Josephson current as in Eq ^j. However, 
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when L > £, the parameters of the normal material be- 
tween the two superconducting leads will affect the tun- 
neling amplitude T kq so strongly that the assumption of 
the constant T kq is not valid any more. In this case we 
write our tunneling Hamiltonian as: 



H = 
H T = 



H R + H L + H n + H 7 
E 

kq 



(A kn Ci <T (r L )Ci k _ <T (r L )d n ,- a (r L )d na (r L ) + 



{A qn Cl{r R )Cl q _ a {r R )d n , 
E n = E(n)d n d n , 



r(r R )d na -(r R ) 



H.C 
(21) 



where dl l (d n ) is the creation(annihilation) operator of the 
electron in the normal material and the index n includes 
all the quantum numbers such as momentum and spin, 
H n is the Hamiltonian of the normal material, E{n) is 
the eigenvalule of H n , A kn is the transmission amplitude 
on the left S/N interface. Ht in Eq. ([21]) describes the 
Andreev reflection and can allow us to consider the effect 
of the phase shifts acquired in the normal region. The 
current in this case is modified to read 



diolr 



(T n g r L (u)T*,g R (u;))n} 



T n = A kn G r (r R ,r L ,n),T:, = A* kn G a (r R ,r L ,n'), (22) 

where G r ^ (r R ,ri,,n, a) is the retarded (advanced) 
Green's function of the electron with energy E{n) and 
spin a in the normal region. When r R — tl = L -C £o> 
G r ( a ) ( rflj r L , n) is close to the Fermi distribution f(E(n)) 
and Eq. (|22|) has the same form to its short junction 
limit Eq. ([2"0"[) . However when r R — rt = L ~ £ , be- 
cause E(n) = —E(n'), G r (r Rl rL 7 n) and G a (r R ,rL,n') 
are not conjugate to each other. As a result, even the 
product G r (r R , rL,n)G a (r R , tl, n') from the trajectories 
with an identical ensemble of scattering centers, which is 
corresponding to the ladder correction, is complex with 
a random phase instead of real. Therefore the ladder 
correction due to the impurity scattering is zero in the 
Josephson junction current. This is very different to the 
calculation of the conductance which is proportional to 
\G r ^(E — Ej)\ 2 where the ladder correction due to the 
impurities is nonzero. This is why impurity scattering 
will not affect the Josephson current. This can be con- 
firmed from the experimental data of the Fraunhofer- 
like interference pattern of the critical current in the 
Josephson junctions^ In those experiments, an oscil- 
lation with an amplitude of almost 100% of the total 
current were observed. 42 Similarly in the ring Joseph- 
son junction with SOI, an oscillation due to the AC 
phase with a larger amplitude than that in the conduc- 
tance experiments is expected. 2 Especially if a one chan- 
nel ring limit is achieved, a 100% oscillation due to the 
AC phase should be observed. When the normal mate- 
rial is semiconducting, a tunable Josephson current con- 
trolled by a voltage can be obtained creating a Joseph- 
son field effect transistor ; 22 i 23 i 24 where the ring Joseph- 
son junction can be switched on and off by tuning the 



AC-phase and carrier concentration. Recently, a voltage 
of —70 V has been reported to switch off a Al/InAs/Al 
Josephson junction^ Our result provides a new possible 
way to control the Josephson current by a gate voltage 
which can be as small as several meV— and much smaller 
than the gate voltage in normal Josephson field effect 
transistors 1 22 i 23 1 2i 



IV. EXPERIMENT PROPOSAL 

In this section, we discuss the feasibility of experi- 
ments to observe AC phase through the oscillation of 
the Josephson current. Our experimental proposal is 
based on the radio-frequency (rf) metho d 36 ! 44 ' 45 which 
is a reliable method to measure the current-phase rela- 
tion (C<i>R) of the Josephson curren t 46 i 47 , especially for 
the small critical current even less than 50 nA4£. There- 
fore this technique makes it possible to detect C$R in a 
S/Sm/S junction although there is a large Schottky bar- 
rier on the S/N interface. Our S/Sm/S junction, depicted 
in Fig. [TJ is incorporated into a SQUID washer of very 
small and well defined inductance L, which is coupled in- 
ductively to a high quality tank circuit in resonance with 
quality Q and coupled strength k Fig. [51 The Joseph- 




Washers 



Tank circuit 

FIG. 8: Proposed experimental set-up. 

son current I s in the ring and the phase difference a of 
the DC current and voltage in the tank circuit have the 
relation 4 ^ 



tan(a) 



fc 2 Q/3/'(0) 



1- 

Ir 



k 2 Ql3 



0f(<t>) ' 

tan a(4>)d(f>, 



(23) 



where (3 = 2irLI c /<5>o is the normalized critical current, 
f(4>) — I s (4>)/I c is the normalized current-phase relation 
and I c is the critical current. Therefore, we can observe 
the oscillation of the Josephson current induced by AC 
phase by detecting the phase difference of the current and 
voltage in the tank circuit. According to the amplitude 
of the Josephson currents and the size of the Josephson 
junctions in Refs. 47.48 .49 we can expect a Josephson 
current up to I c = 1 /zA in the ring shape structure 
whose radius is 1 fim and width of each arm is 300 nm. 
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If we choose the inductance of the superconducting loop 
L = 80pH,— k 2 Q = 0.1, the normalized critical current 
(3 = 2ttLI c /$>o = 0.25, the phase difference a of the 
current and voltage in the tank circuit will be in the 
range (-0.02rad, 0.02rad). 



V. CONCLUSIONS 

In this work we have studied the interplay between the 
critical Josephson current in a ring junction and the AC 
effect. We have calculated the Andreev levels in a Rashba 
type SOI ring system attached to two superconducting 
leads both analytically and numerically. Numerically, us- 
ing the boundary condition of the multiple terminal junc- 
tion in the tight-binding model, we calculate the Andreev 
levels and the Josephson currents in both short and long 
channel length limit. After considering the backscatter- 
ing in the SOI ring, large Schottky barrier on the S/N 
interface, Fermi wavelength mismatch between the super- 
conducting leads and semiconducting ring, and the effect 
of the impurities in the ring, we find that the oscillations 
of the Josephson current due to the AC phase are robust. 
The amplitude oscillation of the Josephson current due to 
the AC effect in the D-C Josephson junction is expected 
to be larger than the conductance oscialltions observed 
in normal ring structures without superconducting ele- 
ments. These results suggest an alternative and likely 
better way to observe the AC phase. Also, since one pe- 
riod oscillation of the AC phase only needs several meV, 
we provide a possible way to create Josephson junction 
FETs controlled by a voltage of the order of meV. 
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APPENDIX A: THE QUASIPARTICLE STATE 
WITHIN THE GAP 



We presents the detail discussion of the quasiparticle 
states in the gap in this appendix. As we know that 
the quasiparticle states in the superconductor can be de- 
scribed by the Bogoliubov equation 



H(r) 
A* 



A(r) 
-H(r) 



u k {r) 
Ufc(r) 



= E 



Uk(r) 
Vk(r) 



(Al) 



where H(r) is the one-electron Hamiltonian, defined as 

h 2 



H(r) 



2m 



■V 2 + U{r) - M , 



(A2) 



where fj, is the electrochemical potential, U(r) is the 
scalar potential and m* is the effective mass. For a ho- 
mogeneous superconductor with A(r)=A and U(r) = 0, 
the solution of Eq lAll can be written as 



ttfc(r) 
Vfe(r) 



exp(i/cr) ( ^° I ,u = cos(|), 



u =sin(-)e-^ ) E 



h 2 h 2 

— — m) 2 + A5A (A3) 



2m* 



where cos 9 = {j^t — [i)/E, <j> is the phase of the su- 
perconductor. |fc| > kf is corresponding to electron- 
like quasiparticle because in this case, cos 9 > 0,9/2 < 
7r/2, cos(#/2) > sin(0/2). For the similar reason, \k < kf\ 
is corresponding to the hole-like quasiparticle. Here only 
the positive energy E is considered. In 1-d case, given 
energy E > A, there are four quasiparticles, two of them 
are the right moving particles and the other two are the 
left moving particles. When considering E < A, accord- 
ing to Eq lA7[ the wave vector k satisfies 



h 2 k 2 
2m* 



[i — ±ie, 



(A4) 



where e = \/A 2 — E 2 > 0. Therefor the wave vector k 
has to take the form 

k — k^ ~\~ iki^ 

where the real part and imaginary part of k satisfy 

fii kiy k^ 

= ±e, 



m 

h 2 

2m* 



(k 2 - k 2 ) - M = 0. 



(A6) 



Given energy E < A, the wave vectors of the quasipar- 
ticles in this case are corresponding to the four solutions 
of Eq|Ml The solution of Eq.jJT]) will be 



V2 



E±ie 



(A7) 



As a result, \u \ = \vq\ = 

If a normal conductor is coupled to a superconductor, 
a unique reflection with energy less than gap, namely 
Andreev reflection, can be observed. To calculate the 
scattering matrix in the interface, we must know the in- 
put and output particles, which are corresponding to the 
right moving and left moving particles separately, in the 
superconducting lead. It is easy to find them when the 
energy is larger than the superconducting gap. However, 
when the energy is less than the superconducting gap, 
the wave vector has to be analytic continuous to com- 
plex plane Eq |A5l It is hard to say which solution is 
corresponding the left moving or the right moving par- 
ticle. We will figure out this problem by some concepts 
from retarded Green's function. 

Retarded Green's function in the position and energy 
space, say G(x, x 1 , E^), is viewed as two traveling wave 
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function outward from a source term 8(x— x'). The oppo- 
site case is so called Advanced Green's function. One way 
to calculate S matrix is to use retarded Green function 33 
but not Advanced Green's function because we are inter- 
ested in the case that if there is an excitation at some 
point, how the excited particle travels away from the 
point of the excitation. This kind of retarded Green's 
function of a fermion is f°° „ e F , c dk where 5 > 0, 
therefore we have a pole where E^ — E + iS. Physically 
we hope that when x > 0, the integral gives us the right 
moving particle and when x < 0, the integral gives us the 
left moving particle. 

Now let us come back to our question. In our discus- 
sion, the energy E is chosen to be positive. As a result, 
Ek must be in the first quadrant of a energy complex 
plane and just a little bit above the real axis. Actually, 
Eq. (|A5[) is a general form of wave vector, ki is infinitely 
small when energy E > A and finite when energy E < A. 
Considering E > A and according to Eq. (|A5[) . the energy 
of the quasiparticles can be written as 



Eh 



h 2 k 2 

2m* 



-fi) 2 + A*A 



2j (^-r - A*)— - — 



2m* 



, fiP 1 k^ . ~ ... , ft Jc 
b-T-M) 2 + A5A + «(—-/*)■ 



2m* 



2m* 



m 



\8) 



here ki is infinitely small and we neglect the kf term. 
Since Ek is in the first quadrant of the complex plane, 
(IMV — P-) h krki > 0. For the right moving electron-like 

V 2m* ' ' m 2 

quasiparticle ^ re , since (f^r — /i) > and k r > 0, ki 
must be also larger than zero, say the wave vector of 
the right moving electron-like quasiparticle of the first 



quadrant in complex wave vector plane. For the similar 
reason we have the results that the left moving electron- 
like quasiparticle ^>i e is in the third quadrant, the right 
moving hole-like quasiparticle is in the second quad- 
rant and the left moving hole-like quasiparticle ^>ih is in 
the fourth quadrant. 

Although we get this result in the case E > 0, it is 
also valid when E < because for the wave vectors, the 
case of E < A is just the analytical continuous of that of 
E > A. The wave functions of these four quasiparticles 
within the superconducting gap are 



®rh = 



*lh = 



\/2 

\/2_ e -ia e ~i4> 
V2 —ia —id) 



V2 
2 

V2 ia —id, 

V2 
2 

V2gi«g-# 



J(\kr\+i\ki\)x 



-i(|fc r |+i|fci|)x 



J(-\k r \+i\ki\)x 



i(\k r \-i\ki\)x 



(A9) 



where a = arccos(i?/A). 

An interesting conclusion should be noticed that in 
both the left and right superconducting leads, exponen- 
tially decay quasiparticles are always the input particles 
and exponentially increase quasiparticles are always out- 
put particles. The velocity of the electrons and the cur- 
rent carried by these electrons can be calculated through 
the velocity operator and electron current operator 



J 



v =< > = ( u* Q Vq ) e 



-i{k r — ihi)x 



I=<y\i\^> = e(u* v *\ e -i(kr-iki)x 



Hk r _OL, 

e e 

m* 



\ p i(k T +iki)x / U Q \ _ 

-ih-i- ( 



-ih-f 

ox 



Vo 

J(k r +iki)x ( u a 



0. 



(A10) 



(All) 



where e is the charge of an electron. The velocity of these decay quasiparticles are zero but the current carried by 
these are not zero. Although Eq. (|A10|) shows that the current is decay, since the quasiparticle will decay to cooper 
pair— which can carry a supercurrent, the total current due to one quasiparticle is / = 



APPENDIX B: SPIN-ORBIT COUPLING IN A TIGHT BINDING MODEL RING SYSTEM 



There have been many theoretical papers talking about the Rashba interaction in a ring system, such as the 
eigenenergy, wavefunction and so on by the analytical method in the continuous case and the exact transmission 
function through the numerical calculation in the tight binding model. However the eigenenergy and wavefunction in 
the tight binding model is still undiscussed. Although this is easy to be derived, we write down the conclusion briefly 
to make our paper more readable. 

The Hamiltonian of the SOI ring system has been given in E dl5l Now we give this tight binding model Hamiltonian 
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of the electron around the point tp = 0. 



H = 



( 


2*o 





-to 


if fJSip/2 








\ 







2t a 


jf p -iS<fi/2 


-to 












-to 




2t 





-to 


it P - iS V>/ 2 








-to 





2i 


it so e^l 2 


-to 












-to 


— if P - i5t P/ 2 


2t 







\ 








—it p-*<V/ 2 


-to 





2*o 


/ 



(Bl) 



cos(|)e- m5¥ ' 
sin(i) 
cos(^) 
sin(2) e ^"- 1 )^ 
y cos(^)e 



X)pin8ip 



J 



/ cos(^)e- mA > \ 

cos(^) 
-sin(^) 
cos(|)e m ^ 
V -sin^e^"^^ J 



(B2) 



The eigenfunction of this tight binding model Hamiltonian is the same as the eigenfunction of the continuous 
Hamiltonian. Acting the tight binding model Hamiltonian Ea lBll on the wave function Eq |B2l and focusing on ip = 
gives us 



2*o(l — cos(n — l)5<p) 2t so sin(n — ^)Sip 
2t so sin(n — \)5ip 2*o(l — cos(nS<p)) 

2t (l - cos(n - -)S(f) cos(-5cp) ( ™^/2) 



sin( 7 /2) 
cos(7/2) 



—2*o sin(n — \)5ip sin(^Sip) 2t so sin(n — h)5tp 

2* so sin(n— ^)8(p 2*o sin(n — ^)Sipsin(^Sip) 

2* (1 — cos(m5ip)) 2t so sin(m + \)8ip \ ( cos(7/2) 



sin( 7 /2) 
cos(7/2) 



2t so sin(m + \)5tp 2* (1 — cos(m + l)5tp) 

n . /-, / l\r \ A c \ ( cos(7/2) 

2* (1 - cos(m + ~)5<p) cos(-5tp) I _ 



-sin( 7 /2) 



When 7 = arctan( 



to sin((5^/2 



~2*osin(m+ i)S(psm(h5(p) 2i so sin(TO+ \)8ip 
2t so sin(m + ^)8tp 2t sin(m + ^)S<p sin(|<5<^) 



■), we have 



cos(7/2) 
-sin( 7 /2) 



(B3) 



(B4) 



E e , n = 2to(l - Vl + (W*o) 2 cos[(n - l/2)8cp + /?]) 
n± = ± arccos (1 



• 



E e , m = 2* (1 - v /T+(W*oFcos[(m + l/2)S<p - 0\) 
E 



m± = ± arccos (1 



(B5) 



(B6) 



where /? = arccos (- 



cos(Sip/2) ^ ^ are Qorresponding to the counterclockwise rotation electron and n_ and 

1 Vi+(t 3 o/to) 2 

m_ are corresponding to the clockwise rotation electron. 



We next show the equivalent between the continous 
limit and the tight bind model by giving the detail proof 
of our statement that the term w so = (2f3/8tp—l) is equal 
1) in the limit a 



to (vT 



0. According to the 



L' Hopital's rule, the term (2f3/Sip~ 1) satisfies 



lim 



2(3 f(a) 



lim 



o Sip g{a) a^o g'{a) ' 

% / cos(8ip/2) 

f(a) = 2arccos( 

= 2 arccos( 



y/l + {tso/to) 2 

cos(a/2r) 



), 



y/l + (ama/h 2 ) 2 ' 



](a) = Stp = a/r. 



(B7) 
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because lim a ^o P = and lim a ^o<V = 0. If we define 

/ \ cos(a/2r) i 

u(a) = , ■■ K ' == , we nave 



^/l+(ama/h 2 ) 2 ' 

darccos(u) 



where 



du 
da 



du 



da 



\J\ - u 2 da ' 



l-v? = 1 - 



cos 2 (a/2r) 



1 + (a 2 m 2 a 2 /H 4 ) 
1 + a 2 m 2 a 2 /h 4 — cos 2 (a/2r) 



1 + a 2 m 2 a 2 /h 4 



a 2 ( 



h 4 

2„,2 



ma 
If 



1 

4r 



J)' 



(B8) 



(B9) 



l/2rsin(a/2r) cos(a/2r)a m a/fi' 



v/1 + a 2 m 2 a 2 /h 4 (1 + a 2 m 2 a 2 /fr 4 ) 3 / 2 

). (BIO) 



fi 4 



By inserting Eq. (|B9IB10p to Eq.(|B8|, we find that 



lim/'(a) 

a— >0 



/ a 2 m 2 1 



^a 2 m 2 /fi 4 + l/4r 2 h 4 4r 2 

(BH) 



1 a 2 ?n 2 
4^2 + -^4- 



Substituting Eq. ([BTT|) to Eq. ([B7)) , we obtain the form 



lim 

a— >0 OiyJ 



/'(a = 0) _ 2yT74^ 2 + a 2 m 2 a/h 4 
9'(a = 0) ~ 



1/r 



1 + 



4a 2 r 



fi 4 



1 + (B12) 



which is exactly the same as our statement. 
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